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, ^ ' Abstract 

\^ [ This analysis is based on the same ideas and numerical inputs of the recent paper 

by Ciuchini et al. on the subject. Some approximations are applied, which make 
'^ ■ analytical calculations applicable in most of the work, thus avoiding Monte Carlo 

t"-. i integration. The final result is practically identical to the one obtained by the more 

^O ■ detailed numerical analysis. Q 
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1 Introduction 



The "2000 CKM-triangle analysis" by Ciuchini et al. |I|] is in my opinion the most accu- 
rate way to use consistently all pieces of experimental and theoretical information relevant 
to infer, within the framework of the Standard Model (SM), the parameters of the uni- 

^ ■ tarity triangle. The basic idea is that beliefs (i.e. probabilities) on the value of each 

input quantity are propagated into beliefs about the output quantities, namely p and f/, 

k> ! and SM related quantities. The propagation is performed using the rules of logic applied 

^ [ to uncertain events, i.e. probability theory, without arbitrariness or "prescriptions". All 

hypotheses are clearly stated and the influence on the results of reasonable variations of 
models and model parameters have been checked and found irrelevant (see Ref. [|^ for 
details). Such a propagation of uncertainty, done by responsible people on physics quanti- 
ties they feel knowledgeable about should not be confused with the kind of mathematical 
games described in Appendix A of Ref. f^. In other words, when we say that the value 
of a quantity has 50% probability to lie in a certain region, given well defined hypotheses 
and with reasonably large flexibility of them, means that we are equally confident that it 
could be there or elsewhere. This is in contrast to "95% CL. results" which should not 
be interpreted as 95% confidence on a certain statement (see e.g. Section 5.2 of Ref. 0|, 
and the Conclusion of Ref. g] for an independent account). 

The motivation of this note is to illustrate a simplified way to perform the analysis 
of Ref. [^. In fact, the Monte Carlo integration used there could seem a bit mysterious. 
As well known, approximations usually loose accuracy, but often help in gaining aware- 
ness. And this is exactly the spirit of this note. Therefore, I have little to add here to 

^ Based on the series of lectures "Statistical methods for frontier physics" given at the VI LNF Spring 
School "Bruno Touscheck", Frascati, Italy, May 2001. A version of the paper with original Mathematica 
eps figures (higher quality, but too large for the arXive) is available at the author's URL. 
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phenomenology, experimental data and inferential framework of the analysis, and I point 
to Ref. |jl| and references therein. 

The note is organized in the following way. In Section the relations between CKM 
parameters are schematically recalled. In the following sections the reweighting effect on 
the points of the p-f] plane due to all pieces of information, is shown in detail. The final 
result is then compared to that of Ref. |I]]. 

2 Uncertain constraints 

The relations (l)-(4) of Ref. [^ can be written in the following way 
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where a,b, . . . e are the constraint parameters, functions of many quantities; the most 
crucial of those determined experimentally are indicated on the left hand side of the 
equations. Note that the order of C3 and C4 is exchanged with respect to Eqs. (3)- 
(4) of Ref. [0]. This is because, since the present information concerning ArUs is of 
different quality with respect to the other quantities, the constraint C4 needs a more 
careful treatment than the others, and it will be introduced after C1-C3. This is also the 
reason why ArUg appears explicitly in the constraint C4. 

In the ideal case all parameters are perfectly known, and the constraints would give 
curves in the p-f] plane. For example, Ci would give a circle of radius ^/a. In other terms, 
due to this constraint all points of the circumference would be appear to us likely, unless 
there is any other experimental piece of information (or theoretical prejudice) to assign 
a different weight to different points. In the ideal case the p.d.f. describing our beliefs in 
the p and fj values would be 

f{p,r^\a)=6{p' + v'-a), (1) 

with the Dirac delta meaning the limit of a very narrow p.d.f. concentrated on the 
circumference. Similar arguments hold for the parameters of the other constraints. Let 
us call them generically p (the fact that C3 depends on two parameters is conceptually 
irrelevant). 

In the real case p itself is not perfectly known. There are values which are more likely, 
and values which are less likely, classified by the p.d.f. f{p)- This means that, e.g. for 
the parameter a, we deal with an infinite number of circles, each having its weight f{a). 
It follows that the points of the (p, f/) plane get different weights. 

The probability theory teaches us how to evaluate f{p,f]) taking into account all 
possible values of p: 

fip,r])= fip,v\p)fip)(^P- (2) 



Now the problem becomes how to evaluate f{p), knowing that each p depends on many 
input quantities Xi, denoted all together with x and described, in the most general case, by 



a joint p.d.f. /(x). In most cases of practical interest, including the one we are discussing, 
each Xi can be considered independent from the others, and the joint distribution simplifies 
to /(x) ^ Yl fi^i)- Calling p(x) the function which relates x to p, f(j>) can be generally 
obtained as 

f{p) = lfi^)Sip-pi^))d^. (3) 

This integral, as well as the previous integrals, is usually performed by Monte Carlo 
techniques. The calculations can be simplified in the following way. 

• First, we rely on the central limit theorem, assuming f{p) to be Gaussian for all 
parameters (but with no constraint on the shape of f{p,fj)\). 

• Expected value, standard (deviation) uncertainty and correlation are basically ob- 
tained by the usual propagation (see Ref. |^ for a similar application). 

• Non linear effects in the propagation have been taken into account up to second 
order, using formulae of Ref. (see this paper for the practical modeling of uncer- 
tainties and treatment of asymmetric cases). 

• The correlation between c and d has been taken into account building a bi-variate 
fic,b). 

As input quantities. Table 1 of Ref. |1[] is used, combining properly (i.e. quadratically) 



the standard deviations. For example, for a I fB^x Bb^ 1 expressed in MeV we obtain 



25^ -|- (20/v3)^ = 27.5, where 20/v3 stands for the standard deviation of a uniform 
distribution of half width 20 MeV. Note that, contrary to what some authors critical 
about Ref. [||] (and the related papers and presentations to conferences) think, I have the 
impression that my colleagues tend to make slightly conservative assessments of uncer- 
tainties. This feeling that I had a priori discussing with them is somewhat confirmed a 
posteriori by the excellent overlap of the partial inference by each constraint (as it will be 
shown in Fig. P) and self-consistency between input parameters and values coming out of 
the inference obtained without the their contribution (see Ref. |]I|]). 

The following results are obtained in terms of expected values and standard deviations 
("p = E(p)±a(p)"): 

a = 0.145 ±0.031 (4) 

b = 0.73 ±0.19 (5) 

c = 3.4 ±0.9 (6) 

d = 1.23 ±0.33 [p{c,d) = 0.76] (7) 

e = 12.7 ±1. 2 ps"^ (8) 

Only the correlation between c and b is relevant, and all others can be neglected, being 
below the 10% level (even that between b and e, related by Am^, is negligible, being 
only ±7%). The radii of the circles given by Ci and C2 are 0.38 ± 0.04 and 0.85 ± 0.11, 
respectively. These radii have the meaning of the sides of the unitarity triangle opposite 
to P and 7, respectively, provided by Ci and C2 alone. 



3 Partial results from Ci, C2 and Cg 
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Applying our reasoning to the constraint given by 
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The contour plot shown in Fig. |l| for fj > 0. 3-D plots are given in Fig. ^ (in all 
figures "rho" and "eta" stand for p and fj). Note that the plot of Fig. |l| is obtained 
slicing /(p, f] I Ci) into 12 iso-p.d.f. contours, and hence the regions shown there have no 
straightforward probabilistic interpretation, since one should make the integrals of the 
p.d.f. inside the region. We leave it as an exercise for the interested readers. 



3.2 Anid 

From the Am^^ we get 
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Figure 1: Contour plot of the p.d.f. of Fig. || for f^ > 0. Note that contours are simply iso-p.d.f. 
levels obtained by 12 contour lines. In order to assign them a probabilistic meaning, one needs to 
evaluate the p.d.f. integrals inside the contours. 
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Figure 2: Probability density function of p and r] obtained by the constraint given by 
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shown by 3-D plots in Fig. ^ 

We can already combine the information from the two constraints, just remembering a 
condition stated in the second paragraph of section ||: In the ideal case we would consider 
the points of the circle equally likely if there were no other experimental or theoretical 
a priori information which would lead to assign different weights to the different points. 
But this is exactly what each partial inference does. Therefore, in order to combine the 
two constraints we just need to multiply the weights of each point and normalize the 
distribution: 

f{p,v\C,,C2)^f{p,v\C,)fip,n\C2). (13) 

The 3-D plot of the normalized distribution is shown in Fig. |^. We see that with these 
two constraints only there is still sign ambiguity, and, in particular, the distribution is 
specular with respect to the f] = axis. 

3.3 \^k\ 

Since the third constraint depends on two strongly correlated parameters, we need to 
consider a joint bivariate Gaussian distribution: 
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Figure 3: f{p,r]) from the constraints given by 
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Figure 4: Probability density function of p and rj obtained by the constraint given by Anid ■ 




Figure 5: Probability density function of p and rj obtained by the constraint given by {ekI- 

with fic = E[c] = 3.4, ac = 0.9, fid = E[(i] = 1.23, ad = 0.33, and pcd = p{c,d) = 0.76. 
The integral 
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can be still evaluated analytically but we omit here the final formula, just giving the 3-D 
plot of the resulting p.d.f. normalized in the region of interest (Fig. ^. 



4 Partial combinations of results (Ci, Ci and C3) 

With this third reweighting, the resulting p.d.f. looses finally sign ambiguities, and be- 
comes rather narrow, with respect to the initial space of possibilities. The 3-D plot is 
shown in Fig. ^ At this point we can evaluate expected value and standard uncertainty 
of the quantities of interest: 



f{p,v\Cl,C2,C3) 



E[p] = 0.19 
alp) = 0.07 

E[f]] = 0.32 
alf]) = 0.04 
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Figure 6: Probability density function and contour plot /(p, 77 1 CI, C2, C3) obtained by the 

, f) and Am^ (see remarks in text and in caption of Fig. [I| about the 



constraint given by 
interpretation of the contour plot). 
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5 Including the experimental information about Am^ 

The treatment of the experimental information about Arris cannot be done by the usual 
uncertainty propagation, just because the simplifying hypotheses (linearization and Gaus- 
sian model for e/ Arris) do not hold. In fact, the likelihood, which has the role of reweight- 
ing the probability, is open, in the sense described in Section 7 of Ref. 0, i.e. it does not 
go to zero at both ends of the kinematical region {Arris = and ArUs = oo in this case) , 
as shown in the top plot of Fig. |^. The reason is simple: in this kind of measurement 
ArUs is not yet incompatible with oo (no oscillation). Though the open likelihood does 
not allow to renormalize the p.d.f. (unless strong priors forbidding high values are used), 
the likelihood can still be used to reweigh the points in the p-f/ plane (see Refs. [§] and |^ 
for other examples and discussions about the R function). 

It is instructive to see how the reweighting of Arris is turned into the reweighting of 
the square radius of circle centered in p = 1 and fj = given by the constraint C4, i.e. 
r^ = e/ Arris. This is illustrated in the central plot of Fig. ^. We see that the reflection 
of i? — *> for Arris < 14 ps~^ essentially kills values above r = 1, resulting in a strong 
sensitivity bound (in the sense of Ref. ) on the angle 7 of the CKM matrix, forced to 
be below 90°. The bottom plot of Fig. |^ shows the reweighting function in the p-f] plane. 

For small values of r^ the reweighting function is divergent, since the whole region of 
high ArUg is squeezed into a small region of r^. This is no serious problem, since these 
points are already ruled out by the other constraints. The cutoff shown in the central 
and bottom plot of Fig. |^ is due to a cutoff at ArUs = 100 ps~^. Note, moreover, that the 
values of Arrig preferred by the data (around 15-20 ps^^) overlap well with the p-f/ region 
indicated by the other constraints. Note also that even if one goes through the academic 
exercise of removing by hand the peak around 15-20 ps~^, chopping R to 1, the effect on 
the values of r above 1 (and hence of 7 above 90°) does not change. 

6 Global combination and results 

Reweighting /(p, fj \ CI, C2, C3) with R function of Arris we get the final result shown in 
Fig. |. 

Expected values and standard deviations are obtained by numerical integration. The 
result is 

E[p] = 0.213 
alp) = 0.042 
fip,r]\Cl,C2,C3,CA) ^ p(p, r/) = -0.12 (17) 

E[f]] = 0.314 
a(r/) = 0.039 

practically identical to p = 0.224 ±0.038 and f/ = 0.317 ±0.040 given in Ref. |I[] using the 
full Monte Carlo integration (obviously, those who pay attention the tenths of standard 
deviations should use the more accurate result of Ref. []I|). 

It is interesting to show partial and global results as contour lines at e~^/^ = 0.61 
of the maximum of the reweighting functions, equivalent to the A(— InL) = 1/2 or 
Ax^ = 1/2 rule^ (I refer to Ref. |jl| for the relation between "standard" methods based 

^Note that the translation of these rules into "confidence" is not straightforward, unless the log- 
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Figure 7: Top piot; likelihood of Arris rescaled to the region of insensitivity. Central plot: 
reweighting factor of r^ = e/Anis- The peak just below 0.2 is an artifact of the cutoff in Arug 
(see text). Bottom plot: same reweighting factor in the p^ plane (note also here the low r 
cutoff). 
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Figure 8: Probability density function and contour plot f{p, rj \ CI, C2, C3, C4) obtained by the 

f], Anid and Arris (see remarks in the text and in caption of Fig. || 
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Figure 9: Top plot: A(lnL) = 1/2 contours from the constraints Ci, C2, C3 and their combi- 
nation. The bottom plot shows the effect of the experimental information on Anig, limited to 
the region of interest. The probability that p and i) are both in the "ellipses" is about 37%. 
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on x^ minimization and the more detailed inferential scheme illustrated there). The top 
plot of Fig. 1^ shows the contour "roads" given by the first three constraints, together with 
the (almost) ellipse of their combination. The probability that the values of p and fj are 
both in the ellipse is about 37%.[| Instead, the projections of the ellipse on each axis gives 
an interval of about 68% probability in each variable. 

The bottom plot of Fig. ^ shows the effect of the constraint C4. First we notice the 
perfect agreement between the Am^^ and Arris roads, indicating that the values of Arris 
suggested by the data are absolutely consistent with the other constraints within the 
Standard Model. Furthermore, the effect of the Am^ on the "ellipse" of the final inference 
is to reshape the left side, increasing the value of p and decreasing its uncertainty, with 
almost no effect on ^, as also shown by the results (|16D and (p!7|). 

7 Conclusion 

The result of Ref. [|l| has been reproduced, using some simplifying assumptions that allow 
most calculations to be performed analytically. The consistency of the pieces information 
on p and f] coming from the different constraints, as well as that relating the constraints 
to the uncertain parameters of the theory (see Ref. fl]] for details), is almost "too good", 
suggesting that the final uncertainties are most likely somewhat over estimated. 

We all are looking forward to the results of the B factories and the next generation of 
experiments on the CKM matrix parameters, hoping to find evidence of new phenomenol- 
ogy that would give new vitality to particle physics. 

It is a pleasure to thank my friends Achille et al. (Ref. |]1[) which have dragged me 
into this business. 
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